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ABSTRACT 

We present high resolution rotation curves of the local group dwarf irregular galaxy 
NGC 6822 obtained with the Australia Telescope Compact Array. Our best curves 
have an angular resolution of 8" or 20 pc and contain some 250 independent points. 
The stellar and gas components of NGC 6822 cannot explain the shape of the curve, 
except for the very inner regions, and NGC 6822 is consequently very dark matter 
dominated. There is no evidence for the presence of a steep density cusp down to 
scales of ~ 20 pc, contrary to the predictions of Cold Dark Matter. 

Key words: galaxies: individual (NGC 6822) - galaxies: dwarf - galaxies: fundamental 
parameters - galaxies: kinematics and dynamics - Local Group - dark matter 



1 INTRODUCTION 

NGC 6822 is a dwarf irregular l ocal group m ember. Located 



at a distance of 490 ± 40 kpc (Mateo 1998), NGC 6822 is 



the most nearby dwarf irregular apart from the LMC/SMC 
system. Due to the small distance the galaxy appears very 
extended on the sky: its optical angular diameter is over 
a qua rter of a degree; the Hi disk measures close to a de- 
gree (|de Blok fc Walter] |2000|). NGC 6822 i s a member of 



68 !2 



an extended cloud of irregulars (Matec 1998) known as the 
Local Group Cloud" . The galaxy has a total luminosity of 
M B = -15.8 ( Hodge et al. 199 jj) an d a total Hi mass of 
1.3 x 10 8 M Q ( de Blok fc Walter||200n"| ), making it relatively 
gas-rich. It is a metal poor galaxy, w ith an ISM abundance 
of about 0.2 Zq (Skillman et al. 1989) and has a star forma- 
t ion rate of ~ 0.06 AJpyr~ (ba se d on H a and FIR fluxes) 
( |Mateo| |l99£j ; [Israel et al] |l996| ). |Hodge| ( |l980| ) found evi- 
dence for i ncreased star f orma tion between 75 and 100 Myr 
ago, while Gallart et al. (1996) showed that the star forma- 



tion in NGC 6822 increased by a factor of 2 to 6 between 
100 and 200 Myr ago. This is consistent with the mostly 
constant but stochastic recent star formation histories often 
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1999). NGC 6822 can be regarded as a rather 



average and quiescent dwarf irregular galaxy. 
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Cosmological numerical Cold Dark Matter (CDM) simula- 
tions predict a specific and universal shape for the dark mat- 
ter mass density d istributions (and therefore the rotation 
curves) of galaxies ( Dubinski fc Carlberg 1991 L This was in- 



vestigated in detail by |Navarro, Ficnk fc White (1996 



1997 



who found that the mass-density distribution in the inner 
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parts of simulated CDM halos could best be described by 
a r _1 power-law. This steeply increasing density cusp to- 
wards the center naturally translates into a steeply rising 
rotation curve. This kind of curve is however not observed 



in late-type disk galaxies (e.g. Swater; 



1999 



Verheijen 



1997 



. The rotation curves of dwarf and Low Surface Brightness 
(LSB) galaxies are shallow and rise linearly, more consis- 
tent with a dark matter distribution that is domina ted by 



a constant-density core with a size of a few kpc (VIoore 
1991 Id e Blok fe McGaugh||l997j |de Blok et al.|fc001b| ; |Blais 



Ouellet pe, Amram fc Carignan [2001 ). Especially in the case 
of LSB galaxies there has been much d iscussion about the 
reality of the ob served shallow curves (Swaters, Madore & 



Trewhe la 



2000). It was argued that systematic resolution 



effects ( "beam-smearing" ) could hide the steep CDM curves 
and lead to an erroneous conclusion that they would not be 
consistent with CDM. 

Subsequent high-resolution follow-up studies, measur- 
ing the rotation curves of LSB galaxies in Ha at a resolution 
of ~ 0.2 to ~ 1 kpc have n ow confirmed the con c lusions de- 



rived from the early data (ie Blok et ; 



Rubin k de Blok 2000; ie Blok et al 



20011 



2001; 



McGaugh 



He Blok & 



Bosma[ 200l|; |Blais-Ouellcttc, Amram fe Carignanj|200l[ [Bo- 



latto et al. [ pool )" Dwarfs and LSB galaxies seem to be dom- 
inated by a dark matter distribution that is best described 
i n the inner parts by a powerlaw p ~ r a with a — —0.2 ±0.2 
( [ie Blok et~aT||2001b| ). 

With a spatial resolution of only ~ 20 pc, the current 
data set enables us to measure the rotation curve of this 
dwarf galaxy at a resolution which is an order of magnitude 
better than the observations described above. Among other 
things, the high resolution makes it possible to investigate 
the effect beam size has on mass models and test whether the 
inferred distribution of dark matter depends on resolution. 



3 OBSERVATIONS AND DATA REDUCTION 

NGC 6822 was observed with the Australia Telescope Com- 
pact Array for 15 x 12 hours in its 375 (1 x 12h), 750D 
(2x 12h), 1.5A (4 x 12h), 6A and 6D (8x 12h) configurations 
over the period from June 1999 to March 2000. A total of 8 
pointings was observed covering the entire Hi extent of the 
galaxy. We used a bandwidth of 4 MHz with a channel sep- 
aration of 0.8 kms" 1 . Additionally, to obtain zero-spacing 
information, NGC 6822 was observed with the Parkes sin- 
gle dish radio telescope using the multibeam receiver in its 
narrowband mode in December 1998. The correlator config- 
uration was identical to that used for the Compact Array 
observations. 

The data were reduced and mosaicked together using 
the miriad data reduction package. Super-uniform weight- 
ing, reducing side lobes in individual pointings prior to mo- 
saicing, was used. The resulting data cubes were cleaned 
with the miriad MOSSDI task. We combined the single dish 
and synthesis data to correct for the missing zero-spacings 
and produced data cubes at various spatial and velocity res- 
olutions. For the rotation curve analysis presented here we 
use the data cubes with channel separations of 1.6 kms" 1 
(with an effective velocity resolution of 1.9 kms -1 ). These 
cubes gave the best compromise between velocity resolution 
and signal-to-noise. 



We analyzed the 1.6 kms 1 data at five different spatial 
resolutions. The synthesized beam sizes used were 96" x 
349.4", 48" x 174.7", 24" x 86.4", 12" x 42.4" and 8" x 28.3". 
The position angle of the beam was 0° at all resolutions. 
We will refer to the 5 different resolution data sets as B96, 
B48, B24, B12 and B08, respectively. Table 1 gives some 
more information on the data sets. We retained the elliptical 
beamshape to obtain the highest possible resolution along 
the major axis. As the PA of the major axis of NGC 6822 
is ~ 110°, as shown later, this means that the resolution 
along the major axis is to a few percent equal to the minor 
axis beam size. The large size of the galaxy compared to the 
beam rules out any beam smearing effects due to non-major 
axis information entering the beam. 



4 MOMENT MAPS AND VELOCITY FIELD 

All subsequent analysis was done using the GIPSY package. 
The low-resolution cube B96 was clipped at the 2a level, and 
remaining noise peaks were removed by hand. Then for B48 
and subsequent resolutions, we used the 2<r-clipped cube of 
the previous resolution as a mask, after which remaining 
spurious noise peaks were removed by hand. This ensures 
a consistent selection of features at all resolutions. Further- 
more using the lower resolution cube as a mask ensures that 
possible extended low-level structures are retained in the 
map. 

The integrated Hi surface density maps were made in 
the usual manner by adding together all clipped channel 
maps. As a second step we then isolated the high signal-to- 
noise (S/N) regions of the maps as follows. For uniformly 
tapered maps in velocity a to t = VNa c h, where a tot is the 
noise in a pixel in the integrated column density map, N is 
the number of channels contributing to that pixel, and a c h 
is the noise in one channel at that pixel. We constructed 
noise maps for each integrated column density map, and 
used these together with the column density maps to isolate 
those pixels in the column density maps where S/N > 10. 
These high S/N maps were used as masks for the velocity 
fields. Figure [l] shows the integrated column density map 
derived for the B12 data. 

There are two standard ways to produce a velocity field. 
One commonly used procedure is to determine the inten- 
sity weighted mean of the velocities along each profile. It 
is known that this method can produce spurious velocities 
at low resolution and low S /N ratios. Given the high resolu- 
tion and S/N of our data set, this is unlikely to be a problem 
here. An alternative method is to fit (a) Gaussian(s) to each 
profile. The latter method usually gives better results, as 
it is less affected by systematic effects due to skewed pro- 
files etc. It is however computationally more expensive which 
becomes important for data cubes of the size we are deal- 
ing with here. The high quality of the current data puts us 
firmly in the regime where the intensity weighted mean gives 
accurate results. 

To test this we have compared both methods for the B12 
data set (Fig. ^) and find that for our data the differences 
are negligible. A histogram of the differences is very well 
described by a Gaussian with an average of —0.3 kms -1 and 
an RMS of 1.0 kms -1 , i.e., both methods produce identical 
results to within better than a channel spacing. The number 
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of pixels where the residuals reach 5 kms" 1 or higher is only 
~ 0.1 percent of the total number of pixels. 

Figure |^ shows a major axis position velocity diagram 
of the B12 data, taken along a position angle of 120°. Over- 
plotted are the velocities as found in the intensity weighted 
and Gaussian velocity fields. The velocity fields produce 
the same results, without large-scale systematic deviations. 
We will therefore proceed to use intensity weighted velocity 
fields in subsequent analysis 



5 DERIVING ROTATION CURVES 

The rotation curves were produced using the GIPSY task 
ROTCUR. For a well-resolved and high S/N data set such as 
the current one, the tilted ring procedure is by far the best 
way to determine the rotation curve. Other methods which 
are often used, such as adjusting the tilted ring parameters 
by hand using position velocity diagrams as a guide (e.g. 
using the GIPSY command inspector), are more subjective 
and results can depend on guesses on e.g. the magnitude 
and effects of beam smearing. 

A look at Fig. [| shows that the high-resolution velocity 
fields contain an extra-ordinary amount of detail. To get 
a good feel for the large scale structure in NGC 6822 we 
started by deriving the B96 rotation curve, and used the 
results for each resolution as the initial estimate for the next 
higher resolution, thus gradually refining the curve. 

We adopted rings with a width of equal to the minor 
axis of the beam. The beam is elongated roughly along the 
minor axis of NGC 6822 with an effective radius of ~ 1.9 
times the minor axis. A ring width equal to the minor axis is 
a good compromise between beam shape and complete sam- 
pling. Most information contributing to the rotation curve 
comes from near the major axis where the beam is narrow. 

In general the procedure involved making a fit with all 
parameters (systemic velocity, centre position, position an- 
gle [PA], inclination and rotation velocity) free. The central 
position and systemic velocity were then fixed and several 
runs with either inclination or PA or both fixed were made 
to find the best model. 

As the minor axis generally provides little information 
regarding the rotation curve we excluded an angle of 30° 
around the minor axis from the fits. One usually also down- 
weights the data around the minor axis by applying a | cos 0\ 
weighting, where 9 is the angle with respect to the major axis 
in the plane of the galaxy. Other alternatives are to apply 
a uniform weighting or a cos 2 6 weighting. We have experi- 
mented with different weighting schemes as well as varying 
the exclusion angle around the minor axis between 15° and 
45°, but found no difference in the curves produced (this is 
due to the large number of independent points that do con- 
tribute at full weight around the major axis). We therefore 
adopt a uniform weighting scheme and a free angle around 
the minor axis of 30° . 

Once a satisfactory curve using both sides of the ve- 
locity field was produced, we also calculated two curves us- 
ing only the approaching and receding sides of the velocity 
field. Comparing the three curves gives information about 
the symmetry of the system. 



5.1 The rotation curves 

For each resolution we derived the curve using the proce- 
dure described above. The results are shown in Table 2 and 
Fig. ^j. The trends found in PA and inclination for the dif- 
ferent resolutions agree well with each other. Only the B96 
data does not show the PA trend as strongly. The variation 
of inclination with radius is only small; only the highest res- 
olutions show weak evidence of a slight warp in the outer 
parts: choosing a constant inclination instead therefore does 
only marginally affect the curves. The trend in PA is pro- 
nounced and real; deriving curves with a constant PA results 
in velocities and inclinations that are inconsistent with the 
data, especially at the higher resolutions. The small-scale 
kinks in the inclination are not physical and mainly due to 
dispersion effects. 

Beyond R ~ 1000" (2.5 kpc) the difference between the 
approaching and receding sides becomes more pronounced 
with increasing resolution. This radius corresponds to the 
edge of the inner Hi disk. The receding curve shows no sig- 
nificant features there, but the approaching curve shows a 
dip in velocity. This coincides with the interface between th e 
main Hi disk and the NW cloud ( |de Blok fc Walteij |200l[ ). 
However, at larger radii the receding curve converges with 
the approaching one, indicating that the galaxy and cloud 
(NW side) as well as the "tails" (SE side) are embedded 
in one symmetrical halo. The global dynamics of the sys- 
tem appear undisturbed despite the morphology of the Hi 
. The small kink at R ~ 150" (0.38 kpc) is caused by a 
number of high velocity dispersion regions that over a small 
radial range happen to be aligned with the tilted rings. At 
R ~ 500" a similar small kink corresponds with the inner 
edge of the large hole. 

Figure ^| shows an overlay of the rotation curve on var- 
ious position-velocity slices. The rotation curve is a good 
representation of the dynamics of the galaxy as a whole. 

For all resolutions we have compared the velocity fields 
with model velocity fields constructed from the tilted ring 
fits, and find no systematic large scale residuals. Figure n 
compares the B12 velocity fields. There are only a few lo- 
calised regions of slightly higher residuals in the higher res- 
olution models, coinciding with regions of high velocity dis- 
persion. Other resolutions give similar results. 



5.2 Final curves 

The error-bars derived from the tilted ring fits are an in- 
dication of the scatter in velocity around the best-fitting 
velocity at each tilted ring. As such they do not take into 
account large-scale asymmetries, differences between ap- 
proaching and receding sides etc. To define more realistic 
errors, we regard the absolute difference between the ro- 
tation velocities of the approaching and receding sides as 
2a errors. We adopt the maximum of the tilted ring uncer- 
tainties and the asymmetry uncertainties as our error-bars, 
with one further modification. In some cases both sides of 
the galaxy are very symmetrical, and these usually are also 
the radii where the tilted ring uncertainties are very small. 
These combine to give very small errors (< O.lkms- 1 ). We 
have therefore imposed a minimum error of half a channel 
width on the curves. After corrections for inclination this 
translates in a minimum error of 0.8 kms -1 . 
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In principle the Hi rotation curves also need to be cor- 
rected for the pressure gradients in the gas to derive the 
true rotation velocities. This correction for asymmetric drift 
is given by 



2 2 



(9 hip 

dkiR ' dkiR 



+ 



where v c is the true rotation velocity, the observed gas 
rotation velocity, a the velocity dispersion in the gas, and p 
the volume density. The median (modal) velocity dispersion 
in the disk of NGC 6822 is ~ 5.8 (~ 5.7) kms" 1 , as mea- 
sured in the 1.6 kms -1 channel spacing data, and does not 
show any obvious radial trends. Assuming a constant scale 
height, we can derive the correction for asymmetric drift. 
We consider the B08 curves where, due to the large gradi- 
ents in the Hi surface density, the corrections are largest. 
The corrected and uncorrected curves are shown in Fig. ^ 
Both curves are very similar, with the largest corrections 
of about ~ 3 kms -1 occuring between 850" and 1050". It 
should be kept in mind that the correction itself is also un- 
certain and depends on assumptions of e.g. constant scale 
height, which may not hold in the outer parts of the gas 
disk. As the uncertainties due to asymmetric drift in the 
inner parts are negligible, ignoring them will not affect the 
results. Based on this analysis, we decided not to apply this 
correction. 



5.3 Dynamical centres 

The central position as derived above varies slightly from 
resolution to resolution. This is understandable given the 
different beam sizes. Fig. ^| shows the central positions over- 
plotted on a B12 column density map of NGC 6822, where 
the error-bars span the major and minor axis FWHM of the 
respective beams. Also plotted are if-band 2MASS (see be- 
low) and J?-band isophotal centers. Due to the large size, 
low galactic latitude and the irregular morphology the opti- 
cal center in NGC 6822 is difficult to determine. From dif- 
ferent fits performed at various isophotal levels we estimate 
an uncertainty of ~ 0.5'. These positions are offset from the 
dynamical center by about 105" or 0.25 kpc. 

Though the dynamical center is by definition the zero- 
point of the rotation curve, it is instructive to re-derive the 
rotation curve under the assumption that the optical center 
is the true center of the galaxy. This is implicitly assumed 
in many emission line observations of rotation curves where 
the optical center of the galaxy is frequently used to line 
up the slit of the spectrograph. We re-derived the B12 ro- 
tation curve fixing the position of the center to that of the 
optical center. We kept the run of position angle and inclina- 
tion identical to that in the original model. Fig. ^ over-plots 
the two curves. The overall shape of the curve is insensitive 
to the precise position of the center. The small-scale differ- 
ences are of the same order of magnitude as the differences 
between the approaching and receding sides of the original 
curve. One should however not exaggerate the importance of 
this offset: had NGC 6822 been at a more typical distance of 
a few tens of Mpc, then the angular size of the offset would 
only be ~ 1" and would have been unnoticeoj^j. 

1 As an aside, we refer ahead to the discussion on central mass- 



6 MASS COMPONENTS 
6.1 Stellar component 

We used a mosaic of full-resolution 2MASS K 3 images to 
derive the surface brightness profile of NGC 6822. The ori- 
entation parameters and the centers of the tilted ring fits to 
the velocity fields were used to measure the average surface 
brightness at each radius. Prior to the ellipse integration 
bright stars were masked out. The average surface bright- 
ness contribution of remaining foreground stars in the field 
away from NGC 6822 was determined and subtracted from 
the profile. 

The K s surface brightness of NGC 6822 is low and at 
large radii the scatter is substantial. At radii beyond 750" 
the signal of the galaxy can no longer be distinguished from 
the sky background. Fig. [ll] shows the azimuthally aver- 
aged surface brightness profile, corrected for inclination, de- 
termined using the B08 tilted ring parameters. The surface 
brightness was corrected for Galactic foreground extinction, 
though this was only a small correction of 0.088 mag in the 
_K"-band. We assumed NGC 6822 to be optically thin in K, 
and only applied a geometric inclination surface brightness 
correction. 

The profile can be well described by an exponential disk. 
As the outer parts of the profile could be marginally affected 
by small uncertainties in the sky background we have fitted 
an exponential disk to the profile at R < 400" and find an 
exponential scale length of 286" or 0.68 kpc. The central 
surface brightness is fj, l (K s ) = 19.8 mag arcsec -2 . The total 
absolute magnitude we derive assuming an exponential disk 
and extrapolating to infinity is Mk = —17.9 or Lk = 3.8 x 
10 8 L . 

For each resolution we determine the surface brightness 
profiles using the relevant tilted ring parameters determined. 
At R > 240" we replace the surface brightness profiles with 
that of the exponential disk fit as indicated in Fig. IxTl The 
rotation curve of the stellar disk was then computed from 
the extended if - band profile following |Sasertario ( 1995 ) and 
Begeman (1987). The disk was assumed to have a verti- 



cal sech distribution w ith a scale height zq = h/6 ( van 
der Kruit & Searle 1981). The rotation curves of the stellar 
component were sampled at the same radii as the rotation 
curves. We assume Y* is constant with radius. While one 
expects some modest variation in Y* with radius in optical 
bands (de Jong 1996), the colour gradients in dwarf LSB 
galaxies in the K band tend to be small, so this effect is not 
likely to be significant (see also Sect. 3.4). 



6.2 Gas component 

We used the orientation parameters of the tilted ring fits to 
determine the azimuthally average Hi surface density pro- 
files using the total column density maps for each resolution. 



density slopes in Sect. 8.1, where we determine the power-law 

within a radius of 0.8 
-0.04 ±0.09. 



slope a of the mass-density profile p < 
kpc. For the B12 case we there find a slope of a 
For the curve presented here, determined using the optical center 
and also using an outer radius of 0.8 kpc, a similar analysis gives 
a = —0.23 ± 0.24. The uncertainties are larger, but the value is 
consistent with a soft core in NGC 6822. The exact position of 
the dynamical center does not significantly influence the results. 
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The Hi surface density profiles are presented in Fig. [H^ In 
derivin g the corresponding rotation curve with the GIPSY 



where ro is the charac teristic radius of the halo and pp the 
characteristic density ( Navarro ; Frenk fc White 1996; Blais 



task rotmod we scaled the surface density profiles by a fac- 



Ouellette, Amram & Carignan 2001). This mass distribution 



tor of 1.4 to take the contribution of helium and metals into 
account. 

We assume that the gas is distributed in an infinitely 
thin disk. This assumption is not crucial. Though one might 
expect dwarf galaxies to have a more puffed up disk than 
a typical spiral galaxy, the shape and amplitude of the gas 
rotation curve depends only very mildly on the thickness of 
the disk. 

To use the B12 curve as an example, for the infinitely 
thin gas disk case the maximum rotation velocity of the gas 
component is 17 kms -1 . When increasing the thickness of 
the disk, this maximum rotation velocity drops only very 
slowly: for an exponential scale height z — 0.5 kpc, it still 
measures ~ 14.5 kms -1 . This value for z is already ap- 
proaching the value of the radial scale length of the stellar 
distribution and likely already too large. We thus adopt the 
thin disk approximation. 



6.3 Halo models 

The dark halo component differs from the previous two in 
that we are interested in parametrising this component as- 
suming some fiducial model. The choice of this model is 
the crux of most of the DM analyses in the literature, and 
many models exist. These can all be broadly distinguished 
in two groups: halo models with a core, and halo models 
with a cusp. An example of the first category is the pseudo- 
isothermal halo, an example of the latter the CDM NFW 
halo. 



6.3.1 Pseudo-isothermal halo 

The spherical pseudo-isothermal (ISO) halo has a density 
profile 



piso(R) = 



1 + {R/Rc) 2 ' 



(1) 



where pp is the central density of the halo, and Rc the core 
radius of the halo. The corresponding rotation curve is given 
by 



V(R) = \jAnG Po R 2 c 



Rc t / R\ 

arr.ta.ri 

\RcJ 



R 



The asymptotic velocity of the halo, Voo, is given by 



Voo = \/4irGp R; 



(2) 



(3) 



To characterise this halo only two of the three parame- 
ters (po, Rc, Voo) are needed, as equation (3) determines the 
value of the third parameter. 



6.3.2 NFW halo 

The NFW mass density di stribution takes the form 
(|Navarro, Frenk fc White] |l996|) 



Pnfw(R) 



po 



(R/r )[l + (R/r W 



(4) 



gives rise to a halo rotation curve 



V(R) = V200 



ln(l + cx) — cx/(l + cx) 
x[ln(l + c) -c/(l + c)] 



1/2 



(5) 



where x = R/R2pp. It is characterised by a concentration pa- 
rameter c = Rzoo/Rs and a radius Jfeoo- These are directly 
related to R s and pi, but are used instead as they are a con- 
venient way to parametrise the rotation curve. The radius 
-R200 is the radius where t he density contrast exceeds 200 , 
roughly the virial radius (Navarro, Frenk & White 1996). 
The characteristic velocity V200 of the halo is defined in the 
same way as i?200- These parameters are not independent 
and are set by the cosmology. 



6.4 Stellar Mass-to-light ratios and weighting 

One of largest uncertainties in any mass model is the value 
of the stellar mass-to- light ratio Y*. Though broad trends in 
T* have been measured and modelled (e.g. Bottema) 



1997 



Bell & de Jong 200C), the precise value for an individual 



galaxy is not well known, and depends on extinction, star 
formation history, Initial Mass Function, etc. The valu e of 
cannot be constrained using the rotation curve alone ( van 



Albada & Sancisi 



|1986| : |Lakc fc Feinswog| p.989| ) and some 



assumptions must be made. We therefore present disk-halo 
decompositions using four different assumptions. 

Minimum disk. This model assumes that the observed 
rotation curve is due entirely to DM. This gives an upper 
limit on how concentrated the dark mass component can 
actually be. It is not a realistic model per se, because it 
ignores the gas disk which is obviously present. However, 
we present it here for ease of comparison with cosmological 
simulations, a large fraction of which tend to only model the 
dark matter component. 

Minimum disk + gas. The contribution of the atomic 
gas (Hi and He) is taken into account, but Y* is assumed to 
be zero. 

Constant T*. Here we choose a value for Y* appro- 
priate for NGC 6822 based on its stellar content, colours 
and likely star formation history. Unfortunately reliable in- 
tegrated colours for NGC 6822 are rare due to its large an- 
gular size. We can however use indirect arguments to arrive 
at a value for Y*. Furthermore we are interested here only 
in the Tf-band value which is rather insensitive to effects of 
extinction and recent star formation. 



Verheijen (1997) presents K'-band rotation curve fits of 



a large sample of HSB and LSB galaxies in the Ursa Ma- 
jor cluster. We consider the rot ation curve fits derived using 
the Bottema disk prescription (Bottema 1997). This recipe 
states that the maximum rotation velocity of the disk is 63% 
of the maximum rotation velocity observed (and is derived 
from measurements of stellar velocity dispersions). The dis- 
tribution of Y+k ratios derived with the Bottema disk for 
the dwarf and LSB galaxies in the UMa sample peaks at 
T t K = 0.4 ± 0.05. If the colours of NGC 6822 are com- 
parable to those found for other dwarfs and LSB galaxies 
(namely B - V ~ 0.5, B - R ~ 8 and V - I ~ 0.8; [h] 
Blok, van der Hulst fc Bothun 1995) we can use Table 3 in 
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Bell & de Jong (200C) to derive an approximate value for 



T* . For their formation epoch model with burst we find val- 
ues for T-kK between 0.35 and 0.40. Similar values are found 
for the infall models and the hierarchical models. From de- 



tailed population synthesis modelling of NGC 6822, Gallart 
et al. (1996) find that the total mass of stars and stellar 
remnant ever formed in NGC 6822 must be ~ 9 • 10 7 Mq. 
Combined with the K-h&nd luminosity derived earlier, this 
gives a value Y*k = 0.24. 

Taking the above into account we adopt a value of 
T*k = 0.35 as representative of the stellar population in 
NGC 6822. 

Maximum disk. The rotation curve of the stellar com- 
ponent is scaled to the maximum value allowed by the ob- 
served rotation curve, but with the restriction that the DM 
density gradient is required to b e flat or negative at all rad ii 
(thus avoiding a "hollow halo" ) ( ^an Albada fc Sancisi 1986 ). 

Each of the rotation curves was fitted using the GIPSY 
task ROTMAS. The program determines the best-fitting com- 
bination of Rc and Voo (for the pseudo-isothermal halo) or 
c and V200 (for the NFW halo), using a least squares fitting 
routine. We assigned weights to the data points inversely 
proportional to the square of their uncertainty. 



that small differences between the tilted ring models also 
play a role. 

We note that at all resolutions we find that (restrict- 
ing ourselves to the ISO model for a moment) the min+gas 
model usually has the smallest value of \ 2 , an d that for all 
cases Xm+g < Xcon- We have tried to constrain the value 
of T* by making a "best fit" to the rotation curves using 
ISO models: i.e. we also let T 4 be a free parameter in the 
fits. Unfortunately the different resolutions make it difficult 
to constrain Y* in this manner. The B96 data is best fitted 
with a stellar component with T*k = 0.10 ± 0.13, while the 
B48 and B24 curves both demand T* A - = 0.31±0.04. For the 
B08 and B12 data such a "best fit" gives unphysical results 
(negative Y*). A constrained fit with T* > yields T* = 
(i.e. the minimum disk+gas case) as best fit for these reso- 
lutions. Though it is clear that T* is difficult to constrain 
using only rotation curve information, the results do suggest 
that the stellar population is dynamically unimportant. 



8 DISCUSSION 



7 RESULTS 

Figures |l^-[l6| show the rotation curves at all resolutions 
plotted for the four assumptions on T* for both NFW and 
ISO halo models. To rule out the possibility that the outer 
parts of NGC 6822, which might be affected by tidal effects, 
affect our conclusions we also present separate fits to just 
the inner part of the galaxy out to to R — 1000" (2.4 kpc) 
or the edge of the main inner disk. Tables ^ and ^ present 
the fit parameters for all models presented here. 

An inspection of Tables [| and ^ shows that the 
goodness-of-fit is in all cases better for ISO than for NFW. 
The only exceptions are some of the maximum disk fits 
where both ISO and NFW are difficult to fit, as NGC 6822 
is not a maximum disk galaxy. A comparison between the 
X 2 values is shown in Fig. [I?]. 

The maximum disk values of T* generally have an un- 
certainty of ± ~ 0.1. Changing the values within this range 
has no discernible effect on the quality of the fits. One point 
of note is that for the B08 case the maximum disk value 
for Y* is less than that for the constant Y* case. The B08 
constant Y* fit thus results in a hollow halo. 

Many of the NFW fits have fit parameters that do not 
make physical sense. The large majority of the fits prefer 
c < and V200 — * 00, which is another way of saying that 
the fit procedure is trying to fit a V ~ R 1 ^ 2 curve to a 
V ~ R curve by stretching it to infinity, thus trying to take 
out the curvature in the model. Where this happened we 
fixed the c parameter to c = 0.1 with only V200 as the free 
parameter. Such a value is still outside the range predicted 
by cosmological simulations where values of c ~ 10 would 
be expected. There is a general trend for the value of the 
maximum disk T* to decrease with increasing resolution. It 
is not clear what the cause of this effect is. The inner slopes 
of the rotation curve and stellar curve change subtly with 
resolution, but not enough to be the sole cause. It is likely 



8.1 Inner slopes 

With the high-resolution rotation curves we can determine 
the inner slope of the mass-dens ity distribution . The ra tio- 



nale behind this is described in de Blok et al 



(2001b). In 



summary, the various halo models make specific predictions 
about the slope of the dark matter mass density distribu- 
tion in the inner parts. If we approximate this distribution 
with a power-law p ~ r a , then the ISO halo predicts a — 0, 
while the NFW halo predicts a = —1. In |de Blok et al 



fl2001bp 



the inner mass density slopes of a large sample of 
LSB galaxies was derived, and the most representative value 
was a = —0.2 ± 0.2, i.e. closer to pseudo-isothermal than to 
CDM. 

Here we derive the slopes of the mass density distribu- 
tion of NGC 6822 at various resolut i ons us ing the inversion 
method described in de Blok et al. j2001b[ ) where we com- 
pute the mass density profile and fit a power-law to the 
mass-density profile. Fig. |l8| shows two examples. We show 
the B96 and B24 mass-density profiles, with the best fitting 
NFW and ISO (minimum disk) models over-plotted. The 
NFW models have central densities that are discrepant by 
a factor of 10 or more. 

We measure the slope within R — 0.8 kpc. This is an 
arbitrary choice, motived by the fact that we want to use 
at least 3 data points to determine the slope of the B96 
curve. Other choices are obviously possible, and may result 
in slightly different values of the slope, but will not affect the 
conclusions. The resulting values are given in Table [s| For 
the B08 data a meaningful slope could not be derived due 
to the large gradients in the rotation velocity that are found 
at small scales. The values are at all resolutions significantly 
different from the NFW values. Note that there is no sys- 
tematic steepening of the slope with resolution, which one 
would expect if the cusp was hidden by resolution effects. 
With a resolution of only tens of parsecs it is not clear that 
there is still room for a cusp in NGC 6822. 
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8.2 Dark matter, feedback and cosmology 

In the previous sections we have made NFW fits to the ro- 
tation curves disregarding any of the predictions that CDM 



makes or c, V200 and their relation. We will compare the 
observed rotation curve with cosmological predictions, it we 
assume that V200 — Knax — 55 kms -1 , we find for ACDM 
that c ~ 9.5 (Navarro, Frenk 
M200 of this halo is M200 



White 1997). The total mass 



= 5.2 • 10 iu Mq. This predicted 
halo is over-plotted on the observed curve in Fig. jl9|, and 
over-predicts the rotation velocity in the inner parts by a 
significant amount. 

To look at the problem from a different perspective, 
let us assume T*=0.35. This gives a stellar mass M* = 
1.3 x 10 s Mq. The total (atomic) gas mass is 1.5 x 10 8 Mq 
(HI+He) and the total observed baryonic mass is 2.8 x 
10 8 Mp. A universal baryon fraction of ~ 0.09 (see e.g. 
White & Fabian 1995) implies a total mass M200 = 3.1 X 



10 9 M Q . For a NFW halo this implies V200 = 21.5 kms" 1 
and R200 = 28 kpc. These values are not consistent with 
the observed curve, and imply (apart from the inability of 
the mo del to fit the data) that either the barvon fraction 



in NG C 6822 differs significantly from the universal value 
or that large amounts of baryons have been expelled, as we 
will discuss below in more detail. 

Le t uo work out the observed matter fractiono in NGC 
6822. We use the ISO models, and dptprminp the amount of 



dark matter within a radius of 5 kpc, the outer edge of the 
Hi disk. This yields a total dark mass out to 5 kpc of 3.2 x 
10 9 Mq. With an observed baryonic mass of 2.8 x 10 s M Q 
we find that NGC 6822 is heavily dark matter dominated 
with M v i s /Mdark(R < 5kpc) = 0.09. This is comparable to 
values found in LSB galaxies and other dwarf galaxies (de 
Blok & McGaugh 1997). The corresponding baryon fraction 
is fb = 0.080, close to the universal baryon fraction. 

The visible matter in NGC 6822 is thus a minor com- 
ponent of the total galaxy system. This makes it harder to 
explain the observed core-dominated dark matter distribu- 
tion as resulting from a NFW halo modified by feedback. 
The term feedback is often used indiscriminately to indi- 
cate star formation and evolutionary processes that affect 
the mass distribution in a galaxy. In fact, there are two dis- 
tinct forms of feedback. One of them is commonly observed 
in galaxies, also in NGC 6822, and usually shows itself as 
small-scale redistribution of mass in the disk. 

The second form of feedback is more catastrophic, 
mostly theoretical, and was introduced to explain the dis- 
crepancy between observed and theoretical dark matter dis- 
tributions. In essence it invokes large-scale and violent star 
formation, resulting in massive blow-outs that drag dark 
matter out due to gravitational interactions, thus destroying 
the primordial CDM cusp. The exact physics of feedback are 
not understood, and the models usually resorts to an empir- 
ical description that is fine-tuned to e.g. fit the Tully-Fisher 
relation or other observational constraints. 

Nevertheless, let us explore the implications of feedback: 
using the amount of energy produced by supernovae and 
making some standard assumptions one can work out an ex- 
pression for the amount of matt er expelled by the effects o f 
supernovae (see e.g. Eq. (A4) in van den Bosch et al. 2000 ). 
Feedback then requires the baryonic mass to drag along a 



similar amount of dark matter (modulo some efficiency fac- 
tor that is unlikely to be much larger than unity). 

Assuming CDM halos, we find for V2 00 — 55 kms -1 
that M 2 q o = 5.2 x 10 10 Mq. Equation A4 in van den Bosch et 
al. ( [2000| ) then yields that the expelled baryonic mass is 3.6 x 
10 y Mq . Assuming that the efficiency factor is unity (every 
unit mass of baryons drags along one unit mass of dark 
matter) we find that the amount of dark matter relocated is 
about equal to the amount of dark matter currently observed 
within the outermost radius of the observed disk. We are 
thus dealing with a major structural re-organisation of a 
galaxy. 

Even if we assume the maximum possible maximum 
disk T*=0.95 we still only find an observed baryonic mass 
of 5.4 x 10 8 Mq (though not physically motivated, this is 
a hard upper limit on the observed baryonic mass). Assum- 
ing an efficiency factor of unity, it seems NGC 6822 must 
have expelled over ~ 6 times its currently observed bary- 
onic mass. 

There has been some discussion in the past whether or 
not stellar feedback can remove a significant frac tion of the 
mass of agalaxy to th e Intergalactic Medi u m, gncdin & 
as well as Mac Low & Ferrara ( 199S ) showed 



Zhao (200 



that this process is unable to explain the observed cores, and 
is also inconsistent with other observational constraints. Ac- 
cording to nume r ical s imulation of stellar feedback by 



Mac 



Low & Ferrara ( 199S ) the potential mass loss is a func- 



tion of the total galaxy mass; for galaxies with gas masses 
< 10 6 Mq a galaxy can potentially destroy iself ('blow-away' 
scenario) whereas for gas masses between 10 6 and 10 y Mq 
some material may be removed ('blow-out', see their Fig. 
1). According to these models we do not expect signifcant 
mass loss in NGC 6822 (M gas = 1.5- 10 s ) due to violent star 
formation. 

From an observer's perspective, there is relatively little 
evidence that mass loss does indeed occur in dwarf galaxies. 
Although strong outflow signatures are observed in some 
cases in dwarf galaxies there are only few cases where outflow 
material may permanently escape the galactic gravitational 
potential (e.g. |Martin| [l998| ). 

A further problem with expelling a lot of mass are the 
global dynamics of the galaxy: if one were to have a huge 
mass loss of order a few times the total galaxy mass, one 
would certainly expect to find some clear evidence in the 
velocity field. The blown-out gas is such a case unlikely to 
end up as an undisturbed regularly rotating Hi disk. Fur- 
thermore, the stellar record shows no evidence for increases 
in the SFR with large fact ors at any time during the past ~ 9 
Gyr (Gallart et al. 1996). It thus seems somewhat surpris- 
ing that if feedback had occurred, it would have managed 
to hide all traces of a violent past. 

In summary, we conclude that potential feedback of sev- 
eral times the total baryonic mass of NGC 6822 is unlikely. 
The observed rotation curve of NGC 6822, combined with 
the A'-band and Hi data thus lead to a picture inconsistent 
with a combination of cuspy halo and feedback. 



9 CONCLUSIONS 

We obtained high resolution rotation curves for the local 
group dwarf irregular galaxy NGC 6822. By fitting both 
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pseudo-isothermal and NFW halo models, it was found that 
the pseudo- isothermal model fits the data successfully. NGC 
6822 is not maximum disk and is dominated by dark matter. 
The highest resolution rotation curve has some 250 indepen- 
dent points and beam-smearing effects are not an issue. The 
position of the dynamical center is unambiguously deter- 
mined. There is no indication of a dark matter cusp down 
to scales of ~ 20 pc. 

Trying to explain the observed curve with a combina- 
tion of feedback and CDM halos leads to several inconsisten- 
cies. Feedback and CDM seem to imply that NGC 6822 has 
lost about 6 times the currently observed amount of baryons. 
There is however no evidence for such violent processes. We 
therefore conclude that (i) the dark matter distribution in 
NGC 6822 is unlikely to have been affected by feedback, and 
(ii) is best described by a model with a constant density core. 
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Figure 1. B12 total column density map. Contours start at 5 Mq pc 2 and increase in steps of 1 Mq pc 2 . The beam is indicated in 
the lower right corner. 
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Figure 2. The intensity weighted velocity field, derived from the B12 data. The black contours run from —55 kms — 1 in the central 
parts to —100 kms -1 in the outer NW parts, in steps of 5 kms~ 1 . The white contours run from —50 kms -1 in the inner parts to +10 
kms -1 in the outer SE parts. The beam is shown in the lower-right corner. 
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Figure 3. B12 data major axis position- velocity diagram taken along position angle 120°. Dotted contour represents — 4cr, full contours 
start at 4c and increase in steps of 2a. The white line shows the values found in the intensity weighted velocity field. The black line 
shows values derived from the gaussian velocity field. 
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Figure 4. Results of the tilted ring fits for the B96, B48 (top row), B24, B12 (middle row) and B08 (bottom-left) resolutions. In each 
panel we show the the position angle of the major axis (top), the inclination (centre) and circular velocity (bottom) as a function of 
radius for the various resolutions. In every top sub-panel the filled circles represent the curve derived using the whole velocity field. The 
open circles represent the curve for the approaching (SE) side, while the triangles represent the receding (NW) side. Overplotted in the 
top and centre sub-panels are the models used to derive these curves. The points in the top sub-panel show the behaviour of PA as a free 
parameter, with inclination fixed to its model value as shown in the center panel. The points in the center sub-panel show the behaviour 
of inclination with PA fixed to the model values as shown in the top panel. The bottom-right panel shows the rotation curves at each 
resolution overplotted. Full line: B96; dotted line: B48; short-dashed line: B24; long-dashed line: b!2; dash-dotted line: B08. 
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Figure 5. Position velocity slices of the B24 data, with the projected rotation curve overplotted. The full contours arc plotted starting 
at 3<r in steps of 2a. The dotted contours show the —3(7 and —5<r levels. There is residual Milky Way emission present at ~ —10 kms -1 . 
The dotted lines represent the center of the galaxy and the systemic velocity. Overplotted is the final curve as derived from the tilted 
ring fits, corrected to the position angle of the slice. Top left: major axis position- velocity slice taken at a positon angle of 115°. Top 
right: minor axis position- velocity slice. Bottom left: position angle of major axis minus 30°. Bottom right: position angle of major axis 
plus 30° . NGC 6822 is, despite the disturbed morphology, a very symmetrical galaxy. 
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Figure 6. Comparison of the observed velocity fields (left column) and the circularly symmetric model velocity fields derived from the 
tilted ring fits (centre column). The white velocity contours start in the center of the velocity field at —50 kms -1 and increase in steps 
of 10 kms -1 . The black contours start in the centre at —55 kms -1 and decrease in steps of 10 km s -1 towards the NW. The right 
most column shows the residual velocity field (observed— tilted ring). The first white contour represents +5 kms -1 and subsequent white 
contours increase with a step of 5 kms -1 . The first black contour shows the —5 kms -1 and subsequent black contours decrease in steps 
of 5 kms -1 . 
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Figure 7. The B08 rotation curve (solid line) and corrected for asymmetric drift (dotted line) 
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Figure 8. Positions of the dynamical centers of the 5 resolutions (crosses) overplotted on the B12 total column density map. Contours 
start at 5 Mq pc -2 and increase in steps of 1 Mq pc -2 . The triangles show the positions of the K band isophotes (top) and R band 
isophotcs (bottom). The inset more clearly shows the relative positions of the dynamical centers. Errorbars span the FWHM beam. The 
centers all coincide to better than a beam width. 




Figure 9. Comparison of the B12 rotation curves derived using the true dynamical center as central position (grey curves), and using 
the ii"-band optical center as central position (black curves). The thin lines represent the curves from the appraching and receding sides 
separately. 
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Figure 10. Hi surface density profiles, corrected for inclination. The various resolutions are distinghuished as follows: B96: full thick 
black line; B48: dotted thick black line; B24: full grey thick line; B12: dotted grey thick line; B08: full thin black line. 
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Figure 11. Azimuthally averaged if-band surface brightness profile. The profile is corrected for inclination and Galactic foreground 
extinction. Spacing between the points is 8". The line indicates the exponential disk fit used at radii R > 240". 
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Figure 12. B96 mass models. The left column shows the fits made using a pseudo-isothermal (ISO) halo. The right column shows the 
NFW halo fits. The top row shows fits made to the entire rotation curve, while the bottom row shows fits to the inner curve (R < 1000"). 
In each panel the following fits are shows: minimum disk (top- left), minimum disk+gas (top-right), constant T + (bottom- left), and 
maximum disk (bottom-right). Fit parameters are given in the sub-panels. An error of "0.00" indicates that this parameter was fixed 
during the fitting. 




Figure 13. B48 mass models. See Fig. 
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Figure 14. B24 mass models. See Fig. [l2[ 
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Figure 15. B12 mass models. See Fig. 
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Figure 17. Reduced chi2 plot, x 2 (ISO) < x 2 (NFW) at all resolutions and T*. The line is the line of equality. Large symbols represent 
the fits to the complete curves, small symbols fits to the inner curves. Filled circles: B96; open squares: B48; crosses: B24; filled triangles: 
B12; and stars: B08. 
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Figure 18. Mass-density profiles for the B96 (left) and B24 (right) models. The long-dashed lines shows the best-fitting minimum disk 
NFW model, the full line the best-fitting ISO model. The dotted lines represent the power-law fits to the data at R < 0.8 kpc. 
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Figure 19. Minimum disk halo for B12 using parameters as required by CDM simulations assuming that V m ax = V200- 
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Table 5. Inner mass density slopes 



Resolution 


a 


Aa 


B96 


-0.13 


0.14 


B48 


+0.02 


0.12 


B24 


+0.22 


0.35 


B12 


-0.04 


0.09 


B08 







High-resolution rotation curve of NGC 6822. 



Table 1. Data cube properties 





B96 


B48 


B24 


B12 


B08 


beam size (") 


349.4 x 96 


174.7 x 48 


86.4 x 24 


42.4 x 12 


28.3 x 8 


pixel size (") 


32 x 32 


16 x 16 


8x8 


4x4 


2.5 x 2.5 



All data cubes presented here have a channel separation of 1.6 kms 1 , and an 
effective velocity resolution of 1.9 kms -1 . 



Table 2. Adopted ROTCUR parameters 





B96 


B48 


B24 


B12 


B08 


XPOS 


19 h 44 m 57 54 s 


19h44">58.76 s 


19 h 44m 58 71 s 


19' l 44 m 58.04 s 


19 h 44 m 58 91 s 


YPOS 


-14°49 m 13.0 s 


-14°49 m 17.6 s 


-14°49 m 21.7 s 


-14°49 m 18.9 s 


-14°49 m 23.8 s 


VSYS 


—53.3 kms -1 


-54.7 kms -1 


—54.7 kms -1 


—54.4 kms -1 


—54.6 kms -1 


PA 


0-1248: 109° 


0-480: 127-105° 


0-480: 130-105° 


0-456: 140-103° 


0-480: 150-106° 




1344-2592: 109-138° 


480-2400: 105-130° 


480-2136: 105-131° 


456-1992: 103-130° 


480-680: 106° 












688-1992: 106-130° 


INCL 


58.7° 


59.9° 


58.7° 


0-1056: 58° 


0-1096: 57° 










1056-1992: 58-68° 


1104-1992: 57-68° 



Italics indicate a range in radius in arcseconds. The pair of numbers following it indicate the linear range in parameter over 
this radial range. 
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Table 3. NFW Rotation curve fits 



Resolution 


model 


c 


Ac 


V200 


AV200 


T* 


V 2 

"■red 


B96 


min 


3.7 


1.2 


130.2 


41.8 





2.73 




mpg 


2.8 


1.4 


157.0 


76.6 





2.49 




con 


0.1 




810.3 


14.1 


0.35 


2.89 




max 


0.1 




721.9 


20.3 


0.95 


5.97 


B96 inner 


min 


0.1 




878.0 


61.6 





3.81 




mpg 


0.1 




1791.7 


56.7 





3.24 




con 


0.1 




636.9 


57.8 


0.35 


3.36 




max 


0.1 




377.5 


65.7 


0.95 


4.36 


B48 


min 


2.1 


1.3 


217.9 


123.4 





3.15 




mpg 


1.8 


1.3 


225.4 


141.0 





2.56 




con 


0.1 




730.6 


14.1 


0.35 


5.91 




max 


0.1 




696.1 


16.9 


0.48 


8.50 


B48 inner 


min 


0.1 




899.6 


26.4 





5.99 




mpg 


0.1 




826.8 


23.6 





4.82 




con 


0.1 




633.6 


33.0 


0.35 


9.38 




max 


0.1 




563.2 


37.1 


0.48 


11.84 


B24 


min 


0.1 




879.5 


11.0 





3.58 




mpg 


0.1 




814.0 


10.0 





2.99 




con 


0.1 




716.7 


15.3 


0.35 


6.95 




max 


0.1 




698.5 


16.7 


0.42 


8.27 


B24 inner 


min 


0.1 




876.5 


27.2 





7.95 




mpg 


0.1 




795.1 


24.7 





6.56 




con 


0.1 




591.3 


29.4 


0.35 


9.31 




max 


0.1 




551.1 


30.5 


0.42 


10.02 


B12 


min 


2.3 


1.1 


193.1 


92.2 





5.74 




mpg 


1.6 


1.4 


237.1 


174.3 





4.42 




con 


0.1 




630.8 


11.5 


0.35 


6.15 




max 


0.1 




630.8 


11.5 


0.35 


6.15 


B12 inner 


min 


0.1 




889.3 


19.2 





10.61 




mpg 


0.1 




792.0 


17.4 





8.69 




con 


0.1 




592.4 


20.1 


0.35 


11.65 




max 


0.1 




592.4 


20.1 


0.35 


11.65 


B08 


min 


0.1 




879.3 


10.2 





8.13 




mpg 


0.1 




797.0 


9.3 





6.73 




con 


0.1 




674.5 


11.7 





10.69 




max 


0.1 




732.6 


10.3 


0.18 


8.24 


B08 inner 


min 


0.1 




898.7 


20.6 





15.52 




mpg 


0.1 




792.7 


19.0 





13.28 




con 


0.1 




590.7 


22.3 


0.35 


18.11 




max 


0.1 




688.2 


20.6 


0.18 


15.52 



"min": minimum disk; "mpg": minimum disk+gas; "con": constant 
T*; "max": maximum disk. Dots (...) indicate the parameter was fixed 
during fitting. 



High-resolution rotation curve of NGC 6822. 



27 



Table 4. ISO Rotation curve fits 



Resolution 


model 


Rc 




P0 


Apo 


Y* 


V 2 

"■red 


B96 


min 


i 


.63 


0. 


.03 


42.4 


1.1 





0.18 




mpg 


l 


71 





04 


37.1 


1.1 





0.24 




con 


2 


12 





06 


27.1 


0.9 


0.35 


0.28 




max 


3 


.10 





.22 


16.7 


1.0 


0.95 


1.08 


B96 inner 


min 


1 


.50 


0. 


05 


45.0 


1.2 





0.05 




mpg 


1 


.47 





.05 


41.3 


1.2 





0.05 




con 


2 


.10 





.17 


26.4 


1.2 


0.35 


0.10 




max 




oo 






11.2 


0.7 


0.95 


0.55 


B48 


min 


1 


40 





03 


47.9 


1.2 





0.80 




mpg 


1 


40 





03 


44.5 


1.3 





1.01 




con 


2 


.03 





.05 


26.2 


0.7 


0.35 


0.67 




max 


2. 


.37 





.07 


21.7 


0.6 


0.48 


0.87 


B48 inner 


min 


1 


22 





06 


53.5 


2.6 





1.19 




mpg 


1 


.12 





.06 


53.3 


2.9 





1.17 




con 


1 


.89 





.21 


27.2 


2.0 


0.35 


1.55 




max 


2 


70 





,55 


20.3 


1.7 


0.48 


1.95 


B24 


min 


1 


51 





03 


44.6 


1.0 





0.99 




mpg 


1 


.58 





04 


39.5 


1.0 





1.31 




con 


2 


.32 





.06 


23.4 


0.6 


0.35 


1.00 




max 


2 


53 





.07 


21.1 


0.6 


0.42 


1.04 


B24 inner 


min 


1 


.31 





07 


49.8 


2.1 





1.06 




mpg 


1 


.20 


0. 


.06 


48.7 


2.2 





0.96 




con 


2 


.27 





.34 


24.1 


1.7 


0.35 


1.60 




max 


2 


88 





,65 


20.6 


1.6 


0.42 


1.84 


B12 


min 


1 


.24 





02 


53.2 


0.8 





0.78 




mpg 


1 


.23 





02 


48.3 


0.9 





0.92 




con 


1 


.84 


0. 


,04 


26.6 


0.6 


0.35 


1.21 




max 


1 


.84 





04 


26.6 


0.6 


0.35 


1.21 


B12 inner 


min 


1 


19 





,04 


55.0 


1.7 





1.53 




mpg 


1 


06 





,04 


54.4 


1.8 





1.35 




con 


1 


.78 


0. 


,15 


27.2 


1.5 


0.35 


2.50 




max 


1 


.78 





,15 


27.2 


1.5 


0.35 


2.50 


B08 


min 


1 


.29 





02 


52.4 


0.8 





1.11 




mpg 


1 


37 





02 


44.6 


0.7 





1.00 




con 


2 


.01 


0. 


,05 


25.8 


0.6 


0.35 


1.84 




max 


1 


.66 


0. 


03 


33.7 


0.6 


0.18 


1.27 


B08 inner 


min 


1 


.43 





05 


49.0 


1.2 





1.81 




mpg 


1 


.39 





05 


44.3 


1.3 





1.73 




con 


3 


.01 





40 


22.2 


1.0 


0.35 


3.06 




max 


1 


.88 


0. 


,11 


31.5 


1.1 


0.18 


2.21 



"min": minimum disk; "mpg": minimum disk+gas; "con": constant 
T*; "max": maximum disk. 



